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Abstract 

The exclusive photoproduction reaction jp — > T(lS)p has been studied with 
the ZEUS detector in ep collisions at HERA using an integrated luminosity of 
468 pb -1 . The measurement covers the kinematic range 60 < W < 220 GeV and 
Q 2 < 1 GeV 2 , where W is the photon-proton centre-of-mass energy and Q 2 is 
the photon virtuality. The exponential slope, b, of the t dependence of the cross 
section, where t is the squared four-momentum transfer at the proton vertex, has 
been measured, yielding b = 4. 31^3 (stat.) l^g (syst.) GeV -2 . This constitutes 
the first measurement of the t dependence of the 7p — > T(1S) p cross section. 
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1 Introduction 



In exclusive photo-production of heavy vector mesons (VMs), J /if) and T, the masses of 
the charm and the bottom quarks provide a hard scale and the process can be described 
by models based on perturbative QCD (pQCD) [HE]. The interaction may be viewed at 
leading order as shown in Fig. [TJ the photon fluctuates into a qq state of small transverse 
size, which interacts with partons in the proton through a two-gluon colour-singlet state, 
forming a heavy vector meson. Thus the cross section is proportional to the square of 
the gluon density in the proton. A characteristic feature of heavy VM photoproduction 
is the rapid rise of the cross section with the photon-proton centre-of-mass energy, W. 
This can be explained through the increasing gluon density with decreasing fractional 
momentum, x oc 1/W 2 (where the x region accessible in heavy-quark production at 
HERA is 10 -4 < x < 10 -2 ). Numerous studies have shown that the dependence of the 
cross section on W can be parameterised as a oc W 5 OH]. Measurements for the J /if) 
meson [5j[6] yielded 5 ~ 0.7. A higher value of 8 « 1.7 has been predicted for exclusive 
photoproduction of T(15) mesons in leading-order pQCD [7J, consistent with the recent 
ZEUS measurement: 5 = 1.2 ± 0.8 [8]. 

Studies of the exclusive photoproduction of light and heavy vector mesons [5] have shown 
that the t dependence of the differential cross section may be approximated in the region 
of small t (\t\ < 1 GeV 2 ) with a single exponential: da/d\t\ oc exp (— b\t\), where t is the 
four-momentum-transfer squared at the proton vertex. The slope parameter, b, measured 
at ZEUS for exclusive J /if) production at W = 90 GeV is b = 4.15 ±0.05 (stat.) ±° f 8 
(syst.) GeV -2 and exhibits a logarithmic variation: b(W) = bo + 2a' ln(W/Wo) 2 , where 
a' w 0.1 GeV -2 . In an optical model approach for exclusive production of VMs, the 
slope parameter b is related to the radii of the proton, R p , and the vector meson, Rvm, 
according to the approximate formula: b ps (_R 2 + Ry M )/A. The value of b measured 
for J / if) production is approximately equal to that expected from the size of the proton 
(6^4 GeV -2 ), in agreement with calculations based on pQCD [9]. This suggests that 
the size of the J / if) is small compared to that of the proton. A similar picture is expected 
in the case of exclusive T(1S) production [TU|[lT]. 

The present paper reports on the first measurement of b in exclusive T(IS') photoproduc- 
tion, observed in the /i + yU.~ decay channel in the kinematic range 60 < W < 220 GeV, 
and complements the previous results [H [121 [13] on T(IS') photoproduction. The data 
correspond to an integrated luminosity of 468 pb _1 , collected in the period 1996-2007. 
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2 Experimental set-up 



In 1998-2007 (1996-1997), HERA provided electron^ beams of energy E e = 27.5 GeV and 
proton beams of energy E p = 920 (820) GeV, resulting in a centre-of-mass energy of 
y/s = 318 (300) GeV. 

A detailed description of the ZEUS detector can be found elsewhere (14J. A brief outline 
of the components that are most relevant for this analysis is given below. 

In the kinematic range of the analysis, charged particles were tracked in the central 
tracking detector (CTD) [15] and, for the data taken after 2001, also in the microvertex 
detector (MVD) [IB]. These components operated in a magnetic field of 1.43 T provided 
by a thin superconducting solenoid. The CTD consisted of 72 cylindrical drift chamber 
layers, organised in nine superlayers covering the polar-angk^l region 15° < 6 < 164°. The 
MVD provided polar angle coverage from 7° to 150°. The transverse- momentum resolution 
for full-length tracks was a(p T )/p T = 0.0058 p T ® 0.0065 ©0.0014/p T , with p T in GeV, for 
data taken before 2001 and a(p T )/p T = 0.0029 p T © 0.0081 © 0.0012/p T , for data taken 
after 2001. 

The high-resolution uranium-scintillator calorimeter (CAL) [17] consisted of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
was subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion (EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections 
(HAC). The smallest subdivision of the calorimeter is called a cell. The CAL energy res- 
olutions, as measured under test-beam conditions, are a(E)/E = 0.18/ y/E for electrons 
and a(E)/E = 0.35/ \[E for hadrons, with E in GeV. 

The muon system consisted of barrel, rear (B/RMUON) [T8] and forward (FMUON) [T4] 
tracking detectors. The B/RMUON consisted of limited-streamer (LS) tube chambers 
placed behind the BCAL (RCAL), both inside and outside the magnetised iron yoke 
surrounding the CAL. The barrel and rear muon chambers covered polar angles from 34° 
to 135° and from 135° to 171°, respectively. The FMUON consisted of six planes of LS 
tubes and four planes of drift chambers covering the angular region from 5° to 32°. The 
muon system exploited the magnetic field of the iron yoke and, in the forward direction, 
of two iron toroids magnetised to 1.6 T to provide an independent measurement of the 
muon momentum. 

1 Electrons and positrons are both referred to as electrons in this article. 

2 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards the 
centre of HERA. The coordinate origin was located at the nominal interaction point for data collected 
before 2001. After 2001 it was redefined as the centre of the CTD. The polar angle, 8, is measured 
with respect to the proton beam direction. 
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The iron yoke surrounding the CAL was instrumented with proportional drift chambers 
to form the Backing Calorimeter (BAC) [19J. The BAC consisted of 5142 aluminium 
chambers inserted into the gaps between 7.3 cm thick iron plates (10, 9 and 7 layers 
in forward, central (barrel) and rear subdetectors, respectively) serving as calorimeter 
absorber. The chambers were typically 5 m long and had a wire spacing of 1.5 cm. The 
anode wires were covered by 50 cm long cathode pads. The BAC was equipped with energy 
readout and position-sensitive readout for muon tracking. The former was based on 1692 
pad towers (50 x 50 cm 2 ), providing an energy resolution a(E)/E = 1.0/ y/E, where E is 
expressed in GeV. The position information from the wires allowed the reconstruction 
of muon trajectories in two dimensions (XY in barrel and YZ in endcaps) with a spatial 
accuracy of a few mm. 

The luminosity was measured using the Bethe-Heitler reaction ep — > cyp with the lu- 
minosity detector which consisted of independent lead-scintillator calorimeter [20] and 
magnetic spectrometer [21] systems. 

3 Kinematics 

The four-momenta of the incoming and outgoing electron and proton are denoted by 
k, k', P and P', respectively. The exclusive reaction under study 

ep — > eTp — > efi + ' pTp (1) 

is described by the following variables (Fig. (TJ top): 

• s = (k + P) 2 , the centre-of-mass-energy squared of the electron-proton system; 

• Q 2 = —q 2 = —(k—k') 2 , the negative four-momentum squared of the exchanged photon; 

• U = (q ' ' P): the fraction of the electron energy transferred to the hadronic final 
state in the rest frame of the initial-state proton; 

• W 2 — (q + P) 2 = —Q 2 + 2y(k ■ P) + m 2 , the centre-of-mass-energy squared of the 
photon-proton system, where m p is the proton mass; 

• M M + M -, the invariant mass of the fi + [J>~ pair; 

• t = (P — P') 2 , the squared four-momentum transfer at the proton vertex, determined 
from the approximate formula: t ~ — (p+ + p~) 2 — (py + p~ ) 2 , where p^ y are the 
components of the transverse momentum of the decay muons. 

The reaction ep — > eTY, (Fig. dj bottom), where Y denotes a hadronic state originat- 
ing from proton dissociation, constitutes an important background. These events mimic 
exclusive T production when the hadrons from proton dissociation remain undetected. 
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Events used in the analysis were restricted to Q 2 values from the kinematic minimum, 
Qmin = m e2/V(l ~~ v) ~ 10~ 9 GeV 2 , where m e is the electron mass, to a value at which the 
scattered electron starts to be observed in the CAL, Q^ax ~ 1 GeV 2 , with an estimated 
median Q 2 value of 10~ 3 GeV 2 . The photon-proton centre-of-mass energy can then be 
expressed as: 

W 2 ^AE p E e y^2E p (E-p z ), (2) 
where E is the energy and pz is the longitudinal momentum of the pair. 

The approximate formula for t introduces dispersion 3 times smaller then that in the 
experimental resolution of this variable after all event selections; approximation (j2J) has 
a negligible effect in the case of W. 



4 Event selection 

Exclusive events in photoproduction were selected online by requiring at least one 

CTD track associated with a F/B/RMUON deposit or with a signal in the BAC consistent 
with a muon. Owing to the inclusion of muon triggers based on signals in the BAC [221123] . 
the rate of recorded dimuon events increased by 17% for a third of the data as compared 
to the previous T(15) analysis [8]. Offline, events were selected as follows: 

• two oppositely charged tracks forming a vertex and no other tracks present in the 
central tracking system; 

• position of the vertex consistent with an ep interaction; 

• both tracks were required to have hits in at least 5 CTD superlayers, to ensure a good 
momentum resolution; 

• transverse momentum of each track pt > 1.5 GeV; 

• \r] + — r]~\< 1.5, where rj^ is the pseudorapidit}{§ of a given track, to suppress Bethe- 
Heitler background (Section [5]); 

• at least one track identified as a muon in the F/B/RMUON or BAC, whenever available 
in a given event [23] ; if not explicitly identified as a muon, the second track had to be 
associated with a minimum-ionising energy deposit in the CAL; 

• 1 7r — a\ > 0.1, where a is the angle between the momentum vectors of fi + and to 
reject cosmic-ray events; 

• invariant mass M M + M - in the range between 5 and 15 GeV; 
3 Pseudorapidity is defined as r\ = — In (tan |). 
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• the energy of each CAL cluster not associated to any of the final-state muons was 
required to be less than 0.5 GeV, in order to be above the noise level of the CAL. It 
implicitly selected exclusive events with an effective cut Q 2 < 1 GeV 2 ; 

• the sum of the energy in the FCAL modules surrounding the beam hole had to be 
smaller than 1 GeV [24] to suppress the contamination from proton-dissociative events, 
ep — > eTY. According to a Monte Carlo study, this corresponds to an effective cut on 
the mass My of the dissociated system originating from the proton, My < 4 GeV; 

• photon-proton centre-of-mass energy 60 < W < 220 GeV and four-momentum- 
transfer squared |t| < 5 GeV 2 . 

The total number of selected fi + fi~ pairs was 2769. The contamination of this sample 
with cosmic ray muons is less then 1%. 

5 Monte Carlo simulation 

The detector and trigger acceptance and the effects due to detector response were de- 
termined using samples of Monte Carlo (MC) events. Exclusive and proton-dissociative 
vector-meson production were simulated with the DIFFVM 2.0 generator [25J. For proton- 
dissociative events, the simulation was supplemented by the JETSET 7.3 MC package [26] . 
For exclusive vector-meson production, s-channel helicity conservation (SCHC) was as- 
sumed. An exponential dependence, e~ 6 '*', was assumed for the differential cross section 
in t with a slope parameter b = 4.5 GeV^ 2 , consistent with the value obtained for exclu- 
sive J/ip electroproduction [SJE]. The W dependence of the 7p — > Tp cross section was 
parameterised as oc W s , with 5 = 1.2 [8]. Electromagnetic radiative corrections associated 
with the decay muons are of the order of 1 % [27] and were not included in the simulation. 

The non-resonant background, consisting of the exclusive and proton-dissociative Bethe- 
Heitler (BH) dimuon events, was simulated using the GRAPE vl.lk MC program |28j . 
After event selection, the contribution of the proton-dissociative events was 25% of the 
Bethe-Heitler MC sample. 

All MC events were generated in the full kinematic range and processed through the 
simulation of the ZEUS detector based on the GEANT programQ [29] and were analysed 
with the same reconstruction and offline procedures as the data. In addition, correc- 
tions [21] of the muon-detector efficiencies determined from a data set consisting of J/ip 
and Bethe-Heitler exclusive production events were applied. 

4 Version 3.13 for the 1996-2000 and 3.21 for the 2003-2007 periods, respectively. 
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6 Determination of the b slope 



The invariant-mass distribution of pairs after applying the selection criteria is shown 
in Fig. [2J The simulated contributions from the Bethe- Heftier (exclusive and proton disso- 
ciative) process and from the T(1S), T(2S) and T(3S) resonances are also presented^]. As 
in the previous paper [8], the BH distributions were normalised to the data in the range 
[5.0-15.0] GeV excluding the [9.0-11.0] GeV mass window where contributions from the 
T resonances are expected. For the determination of the slope parameter for exclusive 
T(1S) production, only events in the mass window [9.33-9.66] GeV were considered. The 
width of the mass window was chosen in order to avoid excessive smearing of the t vari- 
able and to retain a good signal-to-background ratio. According to MC studies, 71% of 
all reconstructed T(1S) events are expected in this window; the relative contaminations 
of T(25) and T(35) states with respect to T(1S) are 1.3% and 0.1%, respectively. The 
contribution from the T(2S) and T(3S) states was neglected for the extraction of the 
slope parameter b. After scanning no cosmic ray muon candidates were found in the 
signal mass window. 

The value of the slope parameter for exclusive Y(IS') production was determined as fol- 
lows: the sum of simulated distributions of all contributing processes was fitted to the 
observed event yields in the signal mass window [9.33-9.66] GeV in the four t bins shown 
in Fig. |3j A binned Poissonian log-likelihood function, ln(L), was used. The expected 
number of Bethe-Heitler background events was fixed to the value obtained from the fi + fi~ 
spectrum outside the signal region as described earlier. Due to insufficient statistics it was 
not possible to evaluate the contribution of proton-dissociative T(IS') events in the final 
sample with the present data. However, the fraction of such events, / p diss, is expected to 
be similar in all diffractive vector-meson production processes [31]. Therefore, the value 
ipdiss = 0.25 ±0.05, determined for diffractive J/ip production [5J, was used. The values of 
the slope parameter for the exclusive and proton dissociative T(1S) production processes 
differ [32]; in the MC the value for the latter was taken to be fe p diss = 0.65 ±0.1 GeV -2 [5J. 

The fit was performed with two free parameters: the slope b and the number of expected 
T(IS') events in the signal mass window. During the parameter scan, the contribution 
of the exclusive T(1S) production to the t distribution was reweighted at generator level 
to the function b ■ exp (— b\t\). The small statistical uncertainties of the MC sample were 
neglected in the fit. The fit yielded: b = 4.3l 2 ;° (stat.) GeV" 2 and 41 ± 10 T(1S) events 
(44% of the events in this mass window). The fit provides a good description of the data; 
the equivalent \ 2 is 0-61 for 2 degrees of freedom. 



5 The ratio of the number of events iV-ms) : ^Vt(2S) : ^t(3S) was fixed in the MC to 0.73 : 0.19 : 0.08 
according to a CDF measurement [30) of the production of T resonances. 
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7 Systematic uncertainties 



The following sources of systematic uncertainty were considered, where the numbers in 
parenthesis correspond to the uncertainties on b in GeV -2 : 

• /pdiss was varied between 0.2 and 0.3, as determined from J/ij) production [5] (io'25); 

• frpdiss was varied by GeV -2 . In addition to the uncertainty from J/ip production 
quoted earlier, the upper variation was extended to the value 6 p di SS = 1-35 GeV -2 
obtained for this parameter when it was also fitted (+01); 

• the contribution of BH events in the mass window [9.33-9.66] GeV was varied between 
55.3% and 56.7%, according to the statistical uncertainty of its normalisation (lo!io)i 

• the fraction of proton-dissociative to all BH events was varied in the range 0.22 to 
0.28 (±0.30). 

Variation of the parameter 5 between 0.7 and 1.7 and variations of the offline selection cuts 
lead to a negligible contribution to the b uncertainty. The total systematic uncertainty 
was determined by adding the individual contributions in quadrature. 



8 Result and discussion 

The slope parameter b for the exclusive production of T(IS') mesons was measured to be 
b = 4.3^ 2 '3 (stat.) 1q 6 (syst.) GeV -2 . A comparison of all HERA measurements of the 
slope parameter b for exclusive light and heavy vector meson production and for deeply 
virtual Compton scattering (DVCS) is shown in Fig. HI This analysis doubles the range 
covered by previous measurements in terms of Q 2 + My M , where Mvm denotes the mass 
of a vector meson. The measured value is in agreement with an asymptotic behaviour 
of this dependence, reflecting the proton radius. This was already suggested by earlier 
measurements and is consistent with predictions based on pQCD models (b = 3.68 GeV -2 ) 



9 Conclusions 

The exclusive photoproduction reaction 7£> — >■ T(15)p was studied with the ZEUS de- 
tector in ep collisions at HERA using an integrated luminosity of 468 pb -1 collected in 
the period 1996-2007. The analysis covered the kinematic range 60 < W < 220 GeV and 
Q 2 < 1 GeV 2 . The measurement of the exponential slope of the t dependence yielded 
b = 4.3l 2 'g (stat.) log (syst.) GeV -2 . This is the first determination of the b parameter 



7 



for T(1S) production. The result is in agreement with expectations of an asymptotic be- 
haviour of the slope parameter as a function of the effective scale present in the process, 
Q 2 + My M . This measurement extends the value of the scale to ~ 90 GeV 2 , the highest 
achieved to date in the measurement of the t-slope parameter for a vector meson. 
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Figure 1: Diagrams for (top) exclusive and (bottom) proton-dissociative vector- 
meson photoproduction in ep interactions. The variables describing the kinematics 
of both processes are introduced in Section^ 
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Figure 2: Invariant mass distribution of pairs. The dashed line shows 

the simulated Bethe-Heitler (BH) (exclusive and proton dissociative) distribution, 
normalised to the data points in the range [5.0-15.0] GeV excluding the [9.0- 
11.0] GeV mass window. Simulated contributions of the T(IS'), T(25) and T(35) 
resonances are shown as a histogram on the mass axis (dashed- dotted line). The 
solid line shows the sum of all contributions. 
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Figure 3: Measured \t\ distribution (full dots) with error bars denoting statistical 
uncertainties. Fitted distributions for simulated events are shown for the Bethe- 
Heitler (dashed line), exclusive T{1S) (dotted line) and proton dissociative T(1S) 
(dashed- dotted line) processes. The solid line shows the sum of all contributions. 
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Figure 4: Comparison of the HERA measurements of the slope parameter b as 
a function of the scale Q 2 + My M for exclusive T{1S) production (the rightmost 
data point), for other exclusive vector-meson production W^U ^SiMS^ and for deeply 
virtual Compton scattering (DVCS) JSffif^Tj /. 
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